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We first explain the pseudogap of high-temperature superconductivity based on an approach of 
quantum optics. After introducing a damping factor for the lifetime r of quasiparticles, the su- 
perconducting dome is naturally produced, and the pseudogap is the consequence of pairing with 
damped coherence. We derive a new expression of Ginzburg-Landau free energy density, in which 
a six-order term due to decoherence damping effect is included. Without invoking any microscopic 
pairing mechanism, this approach provides a simple universal equation of second-order phase tran- 
sition, which can be reduced to two well-known empirical scaling equations: the superconducting 
dome Presland-Tallon equation, and the normal-state pseudogap crossover temperature T* line. 



PACS numbers: 74.25.Dw, 74.62.-c, 71.10.-w 



Besides the unveiled microscopic mechanism, the re- 
markable mysteries of high-temperature superconductiv- 
ity (HTSC) are the superconducting dome, the HTSC 
transition temperature T c as a function of doping, and 
pseudogap in phase diagram [l], 0, [![• It is believed 
that, a correct explanation for the superconducting dome 
and pseudogap is strongly related to the origin of HTSC. 
However, there is no convincing explanation with detailed 
analytical expressions for the superconducting dome and 
pseudogap region. Its origin and physical interpretation 
are fundamental questions in HTSC [l], 0, El- There are 
two most popular physical models on HTSC: Resonating- 
valence-bond (RVB) picture 0, Q and quantum critical 
scenario [f| 0] • The former origins from the idea of spin 
singlet scenario Q and that the pseudogap is a precursor 
to the superconducting state, which is pairing without 
long-range phase coherence [9j. The latter origins form 
the idea that another state of matter, pseudogap phase, 
competes with superconducting gap for the same Fermi 
surface. 

In this paper, we first apply quantum optics approach 
to explain HTSC, and show that the experimental results 
of superconducting dome together with pseudogap can be 
described by a universal equation of second-order phase 
transition. The superconducting dome and normal-state 
pseudogap crossover temperature in HTSC are attributed 
to the finite lifetime of cooper pairs due to decoherence 
damping, which results in localized cooper pairs with 
limited phase coherence. We also present a modified 
general Ginzburg-Landau free energy density (GLFED), 
which includes a sixth-order term due to this decoherence 
damping effect. This phcnomenological two-state model 
does not invoke any microscopic pairing mechanism, but 
can provide a simple universal equation, which can be 
reduced to two well-known empirical scaling equations: 
the superconducting dome Presland-Tallon equation [l(| , 



and the normal-state pseudogap crossover temperature 
T* line [IH. 

Before applying quantum optics approach, we should 
list some commonly accepted facts we ensure for HTSC: 
(1) Superconductivity is due to the formation of Cooper 
pairs with long-range order (l2l . Il3l ] , which are spin sin- 
glet and have d x 2_ y 2 orbital symmetry. (2) Photoemis- 
sion experiments reveal sharp spectral peaks in the ex- 
citation spectrm [idl ]. indicating the presence of quasi- 
particle states. (3) The essential structure is the Cu02 
planes, and the formation of Cooper pairs takes place in- 
dependently within different multiplayers [l5[. (4) The 
electron-phonon interaction is not the principal mech- 
anism of the formation of Cooper pairs [16j |. (5) The 
Cooper pairs are formed from the time-reversed states 
(called as non-pairing breaking) [l7| . 

First we show the BCS Hamiltonian can be expressed 
in quantum optics approach. The Hamiltonian of the 
whole system can be expressed as H = X)i=o > wnere 

p,<r,i q 

is the Hamiltonian of the free system, the Hamiltonian of 
the interaction between quasiparticles and intermediate 
bosons is 

Hi= M P'* ( c p',-.i c P^ia P '- P + H - c ) > ( 2 ) 

p,p',<T,i 

the Hamiltonian of the repulsive interaction between 
quasiparticle and quasiparticle is 



H '2=Yu V.P c p',T.i c ± P 'a,i c -pa,i c P,T,i> ( 3 ) 

pp'i 

and the coupling Hamiltonian between neighboring Cu02 
layers is 
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C p ) fT,i C P.c,i- 



"p,o\i+l C P.c,i 



(4) 



p.ff.l 



and the random Hamiltonian Ha = Y\ „ . h R „ ; . Ha in- 
eludes all the random interaction, such as the in-elastic 
interaction, the random potential caused by doping, and 
so on. This Hamiltonian could introduce a damping rate 
of quasiparticles. Here, ^2' p , ; denotes the sum of quasi- 
particles that participate in Cooper pairing. cp~ a i creates 
a Fermi- type quasiparticle with momentum state p={p X j 
p y ), spin a in the ith Cu02 layer. In the low-temperature 
superconductors, the Cooper pair formed by two elec- 
trons is caused by the electron-phonon interaction. It is 
not true in HTSC as the facts we list above. However, 
there must be the intermediate bosons being changed 
among the quasiparticles in order to induce the attrac- 
tive potential between two quasiparticles. Here we define 
a q as the annihilation operator of the intermediate boson 
with momentum q=(<7 x , %)■ M p i p is the quasiparticle- 
boson interaction strength, V p i p is the quasiparticle- 
quasiparticle interaction strength and U is the coupling 
strength between two neighboring layers. Hi is the en- 
er gy °f the quasiparticle- quasiparticle repulsive interac- 
tion, while H2 is the energy of quasiparticle-boson in- 
teraction, which could cause Cooper pairing. Here, we 
ignore the effect of H3. If \H2\ > \Hi\, the repulsive in- 
teraction will suppress Cooper pairing, for which no su- 
perconductivity occurs. However, if \H%\ < \H\\, many 
Cooper pairs could be formed, which will leads supercon- 
ductivity. 

In conventional BCS theory, the Hamiltonian of 
the electron-phonon interaction H\ is included in the 
electron-electron interaction Hi resulting an effective at- 
tractive potential. However, in HTSC, since there is 
no consensus on the dispersion of quasiparticle and the 
d-type Fermi surface, it is difficult to follow the con- 
ventional BCS approach to investigate HTSC. In this 
case, we include the H2 in H±, and the system Hamil- 
tonian is given by H = Hq + V + H3 + H4, where 
V = ^2' p / a ; /ip, P ',cr,i is the effective quasiparticle-boson 
interaction, with 



C p',er.i C P< CT - i ^P'-P. cr + H.C 



and the effective boson operator A p /_ Pi5 is 



^p' — p,cr 



1 V,p + 

4AV, p C -P'' ff ' iC - 



(5) 



(6) 



where a denotes the opposite a. Since the most impor- 
tant physics is included in Hq + V + H± and the use 
of H3 is only to ensure the long-range order (there is no 
long-range order in the two-dimensional system) , we only 
consider the interaction of the two-dimensional system, 
firstly. 

The interaction Hamiltonian V is composed by many 
two-state systems, in which /ip.p'.o-,; denotes that the two 
states of quasiparticle, |e) = \p',cr, i) and \g) = |p, <r, i) 
are coupled by a boson with moment q = p' p. The 
number of intermediate bosons with momentum q is 

^p'-p ff-^p'-p.ff )■ Here we should emphasize that since 




FIG. 1: The temperature nT c as a function of parameter a and 
density of quasiparticles N q with zero longitudinal damping 
/3 — and AT = 0. For zero a, temperature T is proportional 
to the density of quasiparticles N q . However, in the non-zero 
a case, the dependence of temperature on N q is arc shape, 
and dome shape for larger a, as shown in the inserted firgurc. 



Vp',P > is the repulsive potential and quasiparticle- 
boson interaction leads an attractive potential, we have 



A + A , - 

^V-p,CT^P'-P.CT 



< 



a p'-p a P'-I 



(7) 



which denotes that there is a bias for the Cooper pairing. 
That is to say only the number of bosons larger than a 
certain value can Cooper pairing occurs. It will be seen 



below, h 



p,p',(T,l 



'e,9 



- hR 



-h R - 



h H + M et g (C, 



-h^, a j can be simply expressed as 
H.c). This model has been 



c g A eg 



well investigated in Quantum Optics [19]. In Hciscnbcrg 
picture, (c^ g (t)c e , g (t)} give the number of quasiparticles 
in states \e,g) as a function of time, and (c g '(t)c e (t)) de- 
notes the coherence between two quasiparticle states. On 
the other hand, in the conventional BCS theory, the life- 
time quasiparticle is infinite (or very long). However, 
in HTSC, the random Hamiltonian gives a finite life- 
time r = 7 _1 to quasiparticle since for arbitrary oper- 
ator Q(t) we have j- t Q{t) = ±[Q(t),h R ] = -jQ(t). The 
damping rate 7 includes all the dccohcrcnt factors that 
could destroy the coherence of Cooper pairing and could 
be expressed as a polynome of the density of quasiparti- 
cles, which participate in Cooper pairing, N q (defined be- 
low), 7 = T,jaj(N q y . Here we only consider the first two 
terms, thus 7 can be written as 7 = 7^ + 7^, where the 
transverse damping rate 7t = t^ 1 , which may be caused 
by spin-exchange collision [2fJ, is proportional to N q , and 
7^ = 1 denotes the longitudinal damping rate, which 
is unrelated to N q , but relate to the superconductivity- 
transition temperature T c . 

Considering the statistic, we define the quantum av- 
erage variables as <j> a p(t) = (c+(t)cp(t)} . In Heisen- 
berg picture, the dynamic equation for operator Q(t) 
is ihf t Q{t) = [Q(i),/ip,PW + h R , + and 



ih£ t Q(t) = [Q(i),/i F 



ihjQ(t). Thus, the dynam- 



3 



ics of simplified system could be expressed as 

-70ee + iM(f>ge ~ «^*0e 9 , (8) 

^l^gg - i^ge + «0*^ e g, (9) 

-J(f> eg - ifl {<p ee - (j) gg ) , (10) 



d rh 
dt^ 

d A 



where f2 = 



M, 



■\fnbo 



is the general Rabi fre- 



A tg 



quency. We have used the fact that A eg 
V n bos ~ n nor since (n bos - n nor ) > 1. n bos is the num- 
ber of bosons with moment q = p' — p, and n nor is 
the equivalent boson number caused by the quasiparticle- 
quasiparticlc repulsive potential. Thus, one could see 
that the number of bosons should be larger than the bias 
number n nor . If there is no damping factor, 7 = 0, the 
perfect Rabi oscillation of between states \e,g) is pre- 
sented, which is given by <t> e g{t) = |sin2|0|i. However, 
if 7 7^ 0, 4>eg(t) should be averaged by the lifetime of 
quasiparticle r as 4> e g = — J dt(j)eg{t)e~ t t T , which leads 
to an expression of the second-order phase transition (l9j 



mi 



2 ni 



(ii) 



We define that the quasiparticle density of two states 



are n e and 



"91 



and their sum 



i g = n q is the den- 



sity of quasiparticles with momentum p' and p. Since 
\4>eg\ denotes the average value of transition probability, 
which is also equal with \4> ge \, n e \4>eg\ gives the number 
of bosons radiated by the upper state |e), and n g \4> ge \ 
gives the number of bosons absorbed by the lower state 
\g). By summing all momentum parts and spin, we have 
that N q = J2 p P ' a n q is the total density of quasiparticles 
that participate in Cooper pairing, N bos = ^ p , a rib os 
denotes the total density of intermediate bosons, and 



P,P',C 



is the effective boscn density caused 



by the repulsive interaction. 

Since the transverse damping rate t^ 1 is proportional 
to the number of quasiparticles [2(i| while tT 1 does not, 



thus -±- 

gr T 



aN q and /3 



-, with g 



a is a 



gr L 1 ■— j ft 

proportional parameter, and will be zero for zero tra- 
verse damping rate. In the equilibrium state, we have 



and rip 



nbo 



Thus, the 



n g — ite — 2 9? ojia^ >'>e\'i'eg\ 

relation between temperature T and the density of quasi- 
particles N q could be expressed by an equation of second- 
order phase transition, 

n(T - T c ) = a 2 N 2 q - Q - 2cV) N q + 77AT. (12) 

Here we have used the fact that N nor — Nf, os — (3 2 = 
r)(T — T c — AT), where rj a proportional parameter to 
make equal dimension on the two sides of equal sign. This 
relation is based on the well-known Gingzburg-Landau 
free energy density (GLFED, see below) [HI]. r)AT is a 
temperature bias. T c is the critical temperature, below 



which (T < T c ) the materials are in the superconduct- 
ing state. As we have said, Nb os is the boson number 
contributing the attractive potential, while N nor is the ef- 
fective boson number contributing the repulsive potential. 
Thus, if Nbos < N nor (T > T c ), the repulsive potential 
will suppress the attractive potential. No superconductiv- 
ity occurs in this case. However, the attractive poten- 
tial will suppress the repulsive potential for Nb os > N nor 
(T < T c ), which will leads superconductivity transition. 
In this case, we could get the condition for materials be- 
ing in the superconducting state from Eq. (TTJ 



< -o?N 2 q + 



1 



2a[3 Nq - V AT < rjT c , (13) 



and the critical condition of the superconductivity tran- 
sition 

V T C = -c?N 2 q + ( 7 - 2afi) N q - V AT. (14) 



Figure 1 shows the temperature T c as a function of N q 
for different parameter a with zero longitudinal damping 
rate (3 = and AT = 0. One could see that, for zero 
transverse damping rate, a = 0, the density of quasipar- 
ticles N q is proportional to the temperature T c . However, 
for the non-zero a case, the dependence between T c and 
N q is arc shape, and finally, dome shape. It is very like 
the pseudogap in the diagram of HTSC, in which T* is 
proportional to the doping, and the dependence between 
T c and doping is a dome shape. The detailed discussion 
will be shown in below. By the way, Fig.l coincides with 
the results of /iSR measurements in Ref. [lU . 

Equation (p~2|) is the steady result for the system. We 
could treat it as the consequence of the minimum of 
GLFED f s , if we define an order parameter ipi, of which 
\tpi\ 2 gives the density of quasiparticles N q . In this case 
GLFED could be expressed as 



f s = /„+(A|^| 2 + -B|^i| 4 + C|Vi| 6 ) 



U 



X (^r+i + c.c) , 



(15) 



where A = r) (-AT - T c + T) /2, B = (i - 2 a/3) /4 and 
C = —a 2 /6, and N is the total number of Cu02 layers. 
Below we will see that in the superconductive phase, AT 
is negative, and 4 S> 2a/3. /„ is the free energy den- 
sity of the normal state. As one could see, unlike the 
usual GLFED [2l]], the sixth-order term of \ip\, which is 
completely caused by the damping factors of quasiparti- 
cles, is included. We think the the coupling interaction 
between two neighboring layers less contributes to the 
free energy. AT is a temperature bias. In the special 
case, a = 0, which denotes the lifetime of quasiparticle 
is infinite, Eq. (Ti"5")) will come back to the usual form of 
GLFED. 

From above GLFED one could discuss the sample face 
energy and the response of HTSC in external magnetic 
field. Here we should denote that: Several phenomenal 
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FIG. 2: Energy of pseudogap (blue points) and superconduct- 
ing gap (red points) for a number of HTSCs as measured as a 
function of hole doping x (obtained from Ref. [H|). Lines are 
the consequences derived from our theory with the parameters 
tT 1 = 152 meV, a = 0.27, tiAT = -0.13, y = 0.24, /3 = 0.065 
and \i — 18. Blue line gives the temperature of pseudogap 
T* (a = 0), while red line denotes the superconductivity- 
transition temperature T c . 



theories [1, H2] based on the Ginzburg-Landau equations 
have been established to explain the pseudogap and su- 
perconductivity dome of HTSC. However, these theories 
have two disadvantages: I. One order parameter cor- 
responds to one phase. Thus what one should find is 
the order parameter corresponding to the superconduc- 
tivity phase, but not those for fermion pairing and the 
Bose condensation in RVB picture or the competing or- 
der in quantum critical scenario. II. Pseudogap and the 
(i-density wave (DDW) have not been confirmed to be a 
phase. There is no reason to introduce an order param- 
eter to either. 

In order to compare with the experiments, we should 
write Eq. (fT4I) relating to the doping x. Since N q must 
be a function of x, it could always be expressed as a 
polynome of x, N q = ^2 n ^ n x n . Here we only consider 
the first two terms, N q = v — fix. Thus, Eq.fTJJ could 
be re-expressed as 



T c = T) 



a 2 [i 2 (y-xy + 



2af3j fi{y - x) - rjAT 
(16) 



where y = vj\i. As one could see, for the "localized" mea- 
surements, which do not include the transverse damp- 
ing factor, a = 0, the temperature is linear propor- 
tional to the doping x. We believe it corresponds to the 
temperature T* of pseudogap 0. On the other hand, 
in fact, quasiparticles have finite lifetime, which leads 
non-zero a. The dependence between temperature and 
doping is arc shape, and dome shape for lager param- 
eter a. In this case, Eq. (fl~6)) simultaneously involves 
the well-known two empirical equations [ToL [ill ] : E pg = 
£™ ax (0.27 - a:)/0.22, E sc = £™ ax [l - 82.6(0.16 - x) 2 } . 
Figure [2] shows the comparison between our theoretical 
result and the experimental datum of Ref. [Til ] . They per- 
fectly match each other. We think that the experimen- 
tal techniques [23| , such as angle- resolved photonmission 
(ARP) [lj] and nuclear magnetic resonance (NMR) [IJ], 
only detect the " localized" character of pairs and do not 
include the effective transverse damping factor. This pre- 
dicts T c will increase with small transverse damping fac- 
tor a. 

From above discussion, one could see that the function 
of doping is only to dilute the density of quasiparticles 
(or Cooper pairs) and consequently to increase the effec- 
tive lifetime of quasiparticle (long-range coherence). In 
the low-doped region, the number of quasiparticles is so 
large that their lifetimes are very short, which leads the 
short-range order. At the optimal doping, the number of 
Cooper pairs and their coherence match very well. In the 
over-doped region, the number of Cooper pairs decrease, 
which reduces the superconductivity T c . 

In conculion, we have applied an approach of quantum 
optics to explain the pseudogap in HTSC. By introducing 
the effective lifetime of quasiparticle, the superconduct- 
ing dome is naturally produced. We also derive a new 
expression of GLFED, which includes the six-order term 
of order parameter and could simultaneously give the two 
well-known empirical formulas (Tol . [ill ] . The main results 
of this letter are the modified general GLFED expres- 
sion Eq. (6) and the universal equation of second-order 
phase transition Eq. (8). Despite the simplicity of this 
approach, these general and universal expressions should 
provide new insights into the origin of HTSC. 
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